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Magnetic diffraction patterns which exhibit more than one discrete set of symmetry related peaks 
have been explained by the simultaneous coexistence of more than one magnetic polarisation and 
propagation vector, k, in the scattering amplitude. In this work we suggest using such multi-k 
magnetic structures as prototypical multiple order parameter systems. The examples studied are 
uranium rocksalts, which exhibit magnetic diffraction patterns, the indexing of which requires three 
orthogonal propagation wavevectors. In the cubic phase these systems exhibit, in addition to the 
three sets of symmetry related magnetic diffraction peaks, a set of additional Bragg peaks which 
have been observed by both neutron and x-ray resonant scattering. An interpretation of the neutron 
diffraction form factor of these peaks is presented, identifying these peaks as signatures of quantum 
coherence between the three orthogonal order parameters. 



I. INTRODUCTION 

Experimental probes of quantum correlations are of fundamental interest in developing an understanding of the 
many-body state. Active areas of research include few-atom Bose condensates and strongly correlated electron systems 
cxhibitingjfor example, joint superconducting-magnetic phases. This paper, a sequel to a series of experimental 
papers |2|, y, |^ , addresses the microscopic expression of quantum correlation in bulk condensed matter through a 
direct calculation of the neutron scattering amplitude for both singlc-k (ki) and triplc-k (k.3) diffraction from the 
non-coplanar magnetic configurations of the cubic rock salt compounds USbo.82Teo.i8 and UAs0.sSe0.2- 

Whilst the ki intensities can be interpreted within the conventional cross section, for the k3 terms there are two 
distinct problems. Firstly, given that the triple-k scattering amplitude lies outside the conventional cross section, 
the existence of the long-range-ordered signal needs an explanation; and secondly, the highly structured neutron 
diffraction form factor requires interpretation. The first point has been addressed in earlier work Q , where it is shown 
that the enigmatic triplc-k peaks may be understood as two-site auto-corrclations of a quantum correlator. This 
paper focusscs on inferences that may be drawn from numerical simulations of the relative Bragg neutron diffraction 
intensities. 

The wealth of information contained in the highly anisotropic k3 form factor is interpreted within a generalised 
electromagnetic scattering amplitude founded in Clifford geometric algebra 0, 0. The inferred electronic density 
distribution is an important clue, and constraint, on the physical origin of such k3 peaks Q. 

Assuming that the magnetic interaction is dominated by orbitals of uranium 5/ symmetry consistent with a cubic 
environment, it is found that an interpretation of both the single- and triplc-k form factors is (uniquely) possible via 
the degenerate set of 5/(<5) orbitals of the form 5/(<5) = {(5a; 3 — 3xr 2 )/r 3 , (5y 3 — 3yr 2 )/r 3 , (5z 3 — 3zr 2 )/r 3 }. 

II. EXPERIMENTAL MOTIVATION 

The motivation for this study has arisen from the experimental investigation of non-coplanar, i.e. multi-k, magnetic 
configurations. These have been examined by x-ray resonant scattering and neutron diffraction techniques in the 
uranium rock-salts UAso.sSeo.2, USbo.ssTeo.12 and USbo.s5Teo.15 Q, 0, S d as prototype multiple order parameter 
systems with a view towards understanding, for example, possible joint magnetic-superconducting phases and novel 
spin, charge and dynamic states invoking high-order quantum correlations. 

Multi-k configurations may be understood as an extension of the familiar singlc-k (ki) states found in ferromag- 
nets and other coplanar magnetic structures (e.g. Neel antiferromagnets and helical magnetic structures) which 
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are described by single particle projectors of the state vector yielding magnetic order parameters of the form, 
•Mk Q = Mk Q exp(i6'k Q ), where a is a directional coordinate, Mk a is the polarisation vector and = k Q • r ; + 6 a 
is the relative phase at the i th position. Such ki magnetic conhgurations are specified by one polarisation, Mk a , and 
one propagation wavevector, k a . 

More intricate magnetic structures may be described by a set of state vector projectors, invoking more than one 
order parameter at a given site (a multi-k configuration). For example, a 3-k configuration, which yields 3 inde- 
pendent sets of magnetic Bragg peaks around a given reciprocal lattice vector r corresponding to three (orthogonal) 
propagation wavevectors, requires a set of orthogonal terms, {A4k a = Mk Q exp(i6>k„), M\af, = exp(i6*ka_), Mfcj= 
Mk 7 exp(i#k 7 )} with the polarisation and propagation vectors defining three independent order parameters |l], 0, 0,9 • 

The rocksalts studied exhibit an antiferromagnetic 3-k configuration. The positions and relative intensities of the 
\ti Bragg peaks are described by a longitudinal (H — I — -) moment configuration along each of the principal cubic 
axes, with ordering wave vectors k = (0 1/2) reciprocal lattice units (rlu). The puzzle in these materials is the 
consistent appearance, in both neutron and x-ray resonant diffraction, of a fourth set of Bragg peaks corresponding 
with a Tin symmetry element (c.f. magnetic dipole) with propagation wavevector k3 where k3 = k x + k y + k z . 
Although the k3 peaks are typically between 3 and 4 orders of magnitude weaker than those arising from the primary 
order parameters at ki positions, the body of experimental data suggests that they are not artefacts or multiple 
scattering resonances |l|, |2|, |3j, |4j . Their long range order, temperature, azimuth and electric dipole [Ei-Ti u ] x-ray 
photon energy dependencies attest to them representing a magnetic dipole term with relative phasing (H — I — ) and 
ordering wavevector k 3 = (1/2 1/2 1/2) rlu. 



III. ANALYSIS OF THE ki AND k 3 NEUTRON DIFFRACTION FORM FACTORS 

The k 3 diffraction resonances, requiring simultaneously the vector sum and direct product of the three underlying 
magnetic order parameters to generate the effective dipole and phase respectively, lie outside the conventional two-site 
(level-1) correlation function as sampled by both neutron and [Ei-Ti„] x-ray diffraction techniques 0,0,0- 

Inspection of the neutron data reveals the systematic absence of neutron diffraction peaks both at Qkj =(0 n/2) 
and Qk 3 = (n/2 n/2 n/2) for n = odd integer. The first observation leads to a magnetic structure having a (001) 

dipole moment parallel to ki with (-H )phasing [H 0- 0- 01 : the second suggests an alternating (H — I — ) longitudinal 

(111) moment along k 3 0-01 The two demands are incompatible and a central result of this paper is that such a 
picture is indeed untenable in light of the anisotropy and wavevector dependence of the ki and k 3 neutron diffraction 
intensities. 

Before it is possible to make an analysis of the intensity distribution (form factor), the origin of the scattering 
amplitude must be understood. For the ki peaks there is no problem since they are well explained by the conventional 
neutron diffraction amplitude. However, the k 3 peaks lack a clear interpretation Q, H, y, B> IE @L To this end we 
present a generalised formalism of the scattering amplitude founded in Clifford geometric algebra 0, 0] . 

Fundamental to the approach, which represents spinor matrix mechanics using a geometric algebra based in 
Euclidean 3-D space, is the translation of the set of spinor matrices, {a x , cr y , <t z }, to the set of {e x ,e y ^e z } Eu- 
clidean axes. On introduction of a geometric product, e,ej = 5ij + s/~^- e ijk^k, & n algebra of 8 distinct elements, 
{1; e,, ej, ek] eje^, e^e,; /}, is generated where the product of two orthogonal vectors is called a bivector (bv) and 
J, the unit trivector, is the geometric product of three orthogonal vectors. 

This reformulation of Pauli matrix algebra provides a basis for the systematic analysis of space-time correlations 
of extended systems 0|- A key concept is the projection of the (mesoscopic) multi-particle state- vector ^ onto a 
set of subspaces. For single quasiparticle states, iVi independent subspaces are identified with a set of N% frames of 
quantisation (lcvcl-1 quantum correlation). For joint 'two particle' or lcvcl-2 correlators there are a set of N2 joint 
subspaces and so forth. 

The spinor, ^(x, t), acts as a time dependent field over x serving to rotate the external reference frame to the 
local orientation of the frame of quantisation. A distinct frame of quantisation, i.e. set {e k }, i s associated with each 
(quasi)particle state. An incoherent set of N\ states is then represented by N\ frames of quantisation of arbitrary 
relative orientation, /3 P ' 9 , between the e 3 : e| axes. This is a representation of the disordered state. The transition 
to long-range order may then be represented by one preferred frame of quantisation where the condensation of (3 p,q 
breaks the rotational symmetry of the paramagnet state. The quantisation axis of this global frame, {e k P1 }, abstracts 
the essential geometry of level-1 long range order and may be taken as the microscopic order parameter enabling great 
conceptual and computational simplification. 

Extension to states of multiple electronic order invokes an additional phasing variable, 7 k J , as a measure of the 
phase between elements of the order parameter sets, {e^ Pl } : {e k Pj }. The k 3 relative intensities provide evidence 
that the fourth set of Bragg peaks arises through the ^ phase coherence of the underlying set of single-k order 
parameters 0,0- 
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IV. NEUTRON SCATTERING 



Within Clifford algebra physical observables are expressed as the geometric product of an appropriate state vector 
with a hierarchal set of level-n quantum correlators, {E n , J n } 0. For a level-1, i.e. single quasiparticle representation 
of one has Ex ~ 1 and J x ~ .A/f^ where the ~ sign indicates behaves as. The magnetic cross-section is then 
separated into an interaction together with the quantum-statistical expectation of the (single-particle) magnetic auto- 
correlation ~ (Ji(r, t) Ji(r', <')). Of particular interest to the present work is the level-3 auto-correlation function 
~ (J 3 (r,t)J 3 (r',t')) where J 3 ~ (M z ka + M z kp + - M^M^M^) 0. 

The neutron-electron dipole-dipole scattering amplitude arises from an effective Zeeman field. The magnitude 
depends on the relative separation and spin configuration of the two particles. Written for an electron sitting in the 
neutron field, 

Hi = [giiBWNW ■ v x {-^r ) (!) 

where a and s are the electron and neutron spin operator respectively and the remaining symbols take on their 
conventional meanings ||. Geometrical restrictions, imposed by the form of the electromagnetic interaction, generate 
powerful selection rules. To exploit these it is convenient to express R/\R\ 3 as V(l/|i?|) and subsequently decompose 
(l/|i£|) in plane wave states (we keep the conventional symbol i = T in the expression of harmonic functions) 
giving, 

Hi = in&i ■ [k x (s x k)] — — . (2) 

This expression is microscopic and symmetric both with respect to wavevector and a and s of the electron and neutron 
respectively and may be interpreted as the Hamiltonian of an electron in the effective Zeeman field of the neutron. 
Mapping this two level system to GA gives the electron-neutron dipole interaction term, 

l^ B B k (^ Pj \ij) <- Li B B k [(^Ie° k P ^Ie° p i) - Ie^^Ief Pj ] (3) 

where the left-hand side is the conventional Pauli-matrix notation and the right-hand side the GA equivalent. Since 
the second element of the right-hand bracket has no scalar part, 

» B B k (i,\a% Pj \i,) <-> n B B k (Ie° Pj i>Ie° Pj ^) 

<-> (fi B IB k ) ■ (ipJi^)bv (4) 

where (tpJi^)bv is the bivector projection of the enclosed geometric product 0,0- The interaction may be thus 
viewed as the bivector projection of the (electron) multivector ipJiip^ on the (neutron's) magnetic field. This formalism 
forms a basis for interpreting the scattering amplitude from many electron states at the level-1 projection of the state 
vector. Explicitly, J\ = Ie 3 ■ 3 is the axial vector dual of the eie2 bivector of the j th order parameter OPj and the 
measured autocorrelation function is ([Ief pi ■ (ifj Ji^ )](r ,t)V e ® P1 ' (4'Ji' t l J )](r',t')) wn ' cn ^ s rea( i as the (r, t) «-> (r', t') 
correlation of the bivector projection of the levcl-1 quantum correlator, i.e. (tjjJit/j^), on the quantisation axis of the 
order parameter subspace labelled 1. 

The dipole cross-section of the two-site J\ correlator permits inferences on collinear and single order parameter 
coplanar configurations to be made. However, higher order coherencies within the state vector, such as those present 
in 3-k systems, are not represented, and the conventional formalism yields no clear way forward. 

The strength of the present formalism is then the immediate generalisation to the level-3 correlator, to give the 
interaction basic to the scattering amplitude, (fi B IBk) ■ {^Jz^)bv We suggest to calculate the scattering amplitude 
for the diffraction peaks in the 3-k state through this J3 projection of the state vector and offer, as empirical validation, 
the success of this formalism in accounting for the relative intensities and intricate anisotropies of the k3 peaks (see 
Figure 2 and discussion below). 

In the presence of periodic translational symmetry the J3 correlator is taken to have the following form [(| , 

1 

4 

where the relative inter-order-parameter phase has been set to zero, the three characteristic wave vectors of long range 
order are k^,, k. y and k z respectively, and Ie^ P ^ is the abbreviation of I° p i 'e° P:> '. 

As argued above, the symmetries of the Ji peaks in the 3-k phases of the uranium rock-salts force the representative 
moments in the scattering amplitude to be longitudinal, i.e. parallel to the appropriate ki vector. The j th subspace 
axis of quantisation, e® P ^ , is therefore taken parallel to its ki vector in keeping with the cubic symmetry, and the 
last term of J3 then automatically gives rise to the stationary (Bragg) phase condition for the fourth set of Bragg 
peaks at k3 = k T + k y + k z . 



J 3 = l(I e OP^-r + Ie OP2 e ik v .r + Ie OP3 e ^. r _ j^OPl^.rj^OPi^.rj^OPa^.^ (5) 
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V. CALCULATION OF ki AND k 3 BRAGG PEAKS 

With this as a basis, the relative intensities and anisotropies of the scattering amplitudes for both ki and k3 
Bragg events have been calculated, forming the relevant matrix elements from the 5/ wavef unctions, which split in 
a cubic environment into a singlet and two triplets labelled /?, 5 and e respectively. In the calculations the point 
group symmetries of the ionic orbitals in a cubic environment have been respected; however, the radial extent of the 
orbitals differs from the free ion values in the metallic state. To assess the spatial extent of the 5/ wavefunctions a 
calculation of the ki Bragg intensities has been made in which the set {jo, ji, j&\ of ionic Bcsscl functions compatible 
with cubic symmetry Q are replaced by a single term, exp(— 0.035|Q|) with |Q| in A -1 . This scheme, which replaces 
(jo) by cxp(— 0.035|Q|) whilst setting the (74) and (jg) coefficients to zero, is used to represent the radial extent of 
the uranium ion in rock-salt compounds Since the (jo) coefficient is independent of the orientation of the 

scattering vector with respect to the state vector, the approximation precludes any information being obtained on the 
angular disposition of the wavefunction. 

The results are given as closed circles in the left-hand panel of Figure 1. This single term approximation to 
the {jo, ji, jo\ set reproduces the wave vector dependency of the measured intensities, and enables one, in a semi- 
empirical manner, to address the role of bonding on the radial extent of the uranium valence shell. The lack of 
numerical agreement for the first two data points is attributed to extinction in this large (700 mg) crystal. The figure 
further demonstrates that the marked anisotropies at |Q| ~ 2.9 3.8 A^ 1 and 4.1 A -1 are correctly reproduced 

by the complex vector scattering amplitude of the level-1 projection of the state vector. The right-hand panel of 
Figure 1 is the transform of the integrated intensities of the left-hand panel, i.e. both experimental (open circles) and 
theoretical (closed circles) , to a form factor together with the extinction corrected (open triangles) form factor of pure 
USb [l0|. The latter has been used to set the integrated intensity on an absolute scale. 

Moving now to the k3 peaks, we note that similar, anomalous, |Q| dependent intensities of the k3 reflections are 
observed in UAs . 8 Se . 2 [J 0, & 

and the conclusions drawn are hence likely to be generic to this class of 3-k 
state. From the relative intensities of ki Bragg peaks in UAso.sSeo.2 and USbo.8sTeo.12, taken at |Q| above 2.5 A -1 , 
where extinction corrections appear minimal (see Figure 1), a scaling of the UAso.sSe .2 data (open triangles) to the 
USbo.88Teo.12 data (open circles) has been made. Using this ki scaling, the k3 points of both compounds have been 
amalgamated in Figure 2. 

Direct calculation shows the minimal [(jo) — > exp(— 0.035|Q|); ji = jg = 0] approximation is insufficient to explain 
the intensity-wavevector distribution of the k3 peaks. Therefore the summed contribution from the radial {jo, ji, jo} 
set as determined by the cubic environment is employed at each angular setting. As with the ki peaks, ignoring the 
role of valence shell bonding results in k3 intensities that have poor agreement with the data since the radial extent 
of the wavefunction is artificially restricted. The same ki replacement of the (jo) term is therefore made to allow for 
the bond electron density redistribution in the radial extent of the isotropic component, and, in the same vein, the 
calculations make use of an, empirical, adjustment of the relative (ji) contribution. 

Calculations of the scattering amplitude arising from the {5}, {e} and {/3} sets of symmetry determined orbitals 
in a cubic environment give a ratio of cross sections: 5/(5):5/(e):5/(/3)::30000:100:l in favour of 5/(5). Furthermore, 
the coherent scattering amplitude from the 5/(5) orbitals is uniquely able to reproduce the observed anisotropy of the 
cross section. Most notably this is the case between the {(5/2 5/2 1/2), (7/2 1/2 1/2)} and {(5/2 5/2 3/2), (7/2 3/2 
1/2)} pairs of reflections at |Q| = 3.526 A -1 and |Q| = 3.900 A -1 respectively. The anomalies are captured by the 
5/(5) orbitals whilst the 5/(/3), 5/(e) or a spherically symmetric linear combination of the orbitals all give the |Q| 
= 3.526 A" 1 anisotropy in the opposing sense. Thus, despite systematic absences at Q = (n/2 n/2 n/2 ) for n odd 
integer and the {xyz/r 3 } symmetry, the /3-singlet appears to play no significant role in the 3-k neutron diffraction 
and conceptual images of a moment along (111) directions with x,y,z projections may be misleading. 

A comparison of the relative cross sections for ki and k3 peaks held within the J3 projector for the 5/(5) orbitals 
have a ratio, k3iki :: 50:1. From the empirical intensity ratio of Bragg peaks, k3:ki :: 1:10 3 (see Figures 1 and 2), one 
infers the J\ contribution to the ki peaks is dominant. This observation underlies the successful use of the standard 
cross-section in previous interpretations. The calculated ratios are collated in Table 1. 



J3 correlator 


Relative intensity at Q = r + k n 


5/(5) triplet (k 3 ) 


1 


5/(e) triplet (k 3 ) 


3- 1(T 3 


5/(/9) singlet (ks) 


3- 10" 5 


5/(5) triplet (ki) 


2 • 10" 2 



Table I: Comparative Bragg peak intensities of the J3 correlator. 
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Figure 1: Left-hand panel: ki scattering intensity compared with a first principles calculation (excepting a vertical scaling 
factor) from the 5/(5) orbitals compatible with cubic symmetry. Right-hand panel: data converted to form factor, additional 
triangular points are taken from the observed form factor of pure USb [Tcjl . The latter have been used to set the integrated 
intensity on an absolute scale in the left hand panel (mb/mole). Key: open circles = USbo.8sTeo.12; open triangles = USb; 
closed circles = calculation. The USbo.8sTeo.12 data are not corrected for extinction. 
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Figure 2: k3 scattering intensity compared with a first principles calculation (excepting a vertical scaling factor) from the 
degenerate 5/(<5) orbitals compatible with cubic symmetry. Key: open circles = USbo.ssTeo.12; open triangles = UAso.sSeo.2; 
closed circles = calculation. The {(5/2 5/2 1/2), (7/2 1/2 1/2)} and {(5/2 5/2 3/2), (7/2 3/2 1/2)} pairs of reflections at |Q| 
= 3.526 and |Q| = 3.900 A -1 respectively have been separated by 2% in their abscissae to illustrate the anisotropy. The data 
are set on an absolute scale as in Figure 1. 
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VI. IMPLICATIONS AND CONCLUSION 

The Clifford level-3 quantum correlators, interpreted as projectors of the mesoscopic state vector, yield a scattering 
amplitude comprising both a sum and geometric product of three independent order parameters [f| , permit a consistent 
interpretation of the observed peaks in terms of three coherent scattering amplitudes coming from the 5/(<5) triplet 
state. The magnctic-dipolc terms of (i/}J n ip^)bv', ft = 1,3 are then extracted by the projection effected by the scattering 
amplitude onto 3 orthogonal subspaces which may be identified by their respective order parameters. The individual 
and summed scattering amplitudes over the 5/(<5) subspaces then determine the ki and k3 cross sections respectively. 
The present calculations of the neutron dipole form factor show that the origin of the 3-k Bragg peaks is consistent 
with selection of the Je* Jejs Je* component of the J3 quantum correlator through the wavevector selectivity of 
coherent diffraction, i.e. fulfilling the Bragg condition at k3 = (kkk). Further, the overwhelming relative amplitude 
of the 5/(5) triplet and its unique ability to explain the observed anisotropy indicate that it is the most plausible 
candidate as the orbital basis to understand the k3 diffraction peaks within this formalism. 

The existence and anisotropy of k 3 Bragg reflections give unique evidence for the space-time long range order of 
a level-3 projection of the state vector. The advantage of the magnetic system studied here is the direct observation 
afforded by neutron diffraction. Associating a set of 5/(5) orbitals with each uranium site, the present work suggests 
that: (i) valence shell bonding gives a radial density redistribution within the cubic point group symmetry, and (ii) 
the electronic system sits in a spatially correlated coherent superposition of the set of {±<5[x], ±5[y], ±<J[z]} orbitals. In 
particular, the k3 null result at (ft/2 ft/2 ft/2) for n odd integer leads to the conclusion that the set {M^ a , ■Mu^ 1 -^k 7 } 
of vector projected order parameters are aligned with the {±Q X , ±Q y , ±Q Z } axes respectively. Thus, the relative phase 
of the order parameter set is ferromagnetic and it follows that {M^. , M.^ } and {±Q X , ±Q y , ±Q Z } share the 

same chirality. 

There is a distinction between the scattering state vector, vp, and the scattering amplitude as the bivector projection 
of (\I>J„^'). For ki peaks, whether they derive from J\ or J3 quantum correlations, the difference is not apparent. 
It comes to light for k3 peaks which derive from spatially extended correlations of the product term associated with 
the level-3 quantum correlator. The J3 correlator projects three scattering amplitudes associated with the three 
independent Pauli subspaces. 

These three amplitudes give rise to two distinct sets of diffraction peaks associated with the two level- 1 and level-3 
projections of the state vector. It is suggested that the two sets may be identified with two types of electronic order, 

(antifcrro)magnetic (-1 — I ) ki and (antiferro)chiral (-1 — I — ) k3 order respectively. The notion of ki as primary and 

k3 secondary order parameters, based on the respective thermal evolution of the associated sets of Bragg reflections, 
may in this respect be a result of the temperature dependencies of the J\ and J3 correlators respectively. 

It would appear difficult for arguments based on semi-classical approximations to level- 1 projections of ^ to capture 
the essence of the triple-k quantum state, and by implication, of the non-coplanar magnetic state in general. The 
relevant correlations appear to lie outside of such a formalism. A similar limitation may confound intrinsically level- 1 
approaches to problems such as the joint superconductor-magnetic phases and quantum critical point scenarios, which 
are of great topical interest. 

The order parameter concept is general. A particular advantage of the present realisation of multiple order param- 
eters as magnetic states, where "J is encoded as a rotor giving the space-time local orientation of a spin-vector, is the 
clear geometric interpretation. Wc hope that this work will stimulate further research in multiple order parameter 
systems. The role of the quantum correlators of the mesoscopic state vector may be of direct interest in the study 
of more complex multiple order parameter systems such as magnetic superconductors and mutual charge, orbital 
and magnetic ordered states. In particular the present study motivates the search for experimental evidence of both 
further states of long range order in level-n quantum correlators and the existence of novel states that may be revealed 
through (J1J3) cross-correlators. 
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